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1  Introduction 

The  modeling  strategies  for  Large  Eddy  Simulation 
(LES)  of  premixed  flames  have  been  based  mainly  on 
G-Eq  (Peters  2000;  Moureau  et  al.  2009),  flame  sur¬ 
face  density  (Hawkes  and  Cant  2001;  Lecocq  et  al. 
2011)  and  laminar  flamelets  combined  with  various 
forms  of  presumed  probability  density  function  (pdf) 
to  express  flame  filtering  (Kim  et  al.  1999;  Domingo 
et  al.  2005;  Fiorina  et  al.  2010;  Moureau  et  al.  2011; 
Lecocq  et  al.  2011;  Nambully  et  al.  2014a;  Nam- 
bully  et  al.  2014b).  Artificial  thickening  of  the  flame 
front  (Colin  et  al.  2000;  Kuenne  et  al.  2012;  Boudier 
et  al.  2008)  is  another  successful  alternative. 

A  fully  new  modeling  strategy  is  discussed  in  this 
paper,  in  which  an  approximate  deconvolution  is  com¬ 
bined  with  explicit  filtering  to  express  the  sub-grid 
scale  (SGS)  terms  of  scalar  fields. 

This  idea  goes  back  to  the  fundamentals  of  LES 
grounded  on  physical  space  filtering,  instead  of  intro¬ 
ducing  an  ad-hoc  treatment  of  the  phenomena  which 
are  unresolved  by  the  grid. 

The  proposed  method,  called  ADEF  (Approximate 
Deconvolution  and  Explicit  flame  Filtering),  may  be 
applied  to  any  form  of  chemical  description,  Arrhenius 
expressions  or  tabulated  chemistry. 

First,  the  method  is  applied  to  a  canonical  lam¬ 
inar  one -dimensional  filtered  flame.  Then,  LES  of 
a  three-dimensional  turbulent  premixed  jet  flame  is 
performed  and  results  are  compared  against  measure¬ 
ments  by  Chen  et  al.  (1996). 


2  ADEF  Modeling 


Approximate  deconvolution 

A  Taylor  expansion  is  written  for  the  scalar  field  in 
every  ‘j’  direction: 
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Repeated  indices  implies  a  sommation  over  the  three 


directions.  A  Gaussian  filter  is  applied  to  relation  (1), 
all  terms  with  odd  power  of  the  position  are  eliminated 
due  to  the  filter  symmetry.  For  turbulent  scalar  signals, 
the  fourth  order  terms  were  not  found  to  play  a  major 
role,  and  the  filtering  may  be  achieved  with  only  the 
second  order  derivatives  (Katopodes  et  al.  2000): 

«fc»)-*te«)  +  fj3^+0(A*)-  (2) 

The  relation  (2)  may  be  used  recursively,  and  rearrang¬ 
ing  leads  to  an  approximate  and  discrete  deconvolution 
operation  : 

A2  d2T5 
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A  density  weighted  scalar  deconvolution  may  also  be 
introduced  as: 
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The  relations  (2)  and  (3)  may  be  used  in  explicit  or 
implicit  numerical  formulations. 

Application  to  flames 

Filtering  the  (^-balance  equation  of  a  scalar  used  to 
describe  a  flame  front  leads  to: 

^  +  V  •  (palp)  =  V  •  (pZ4@V^)  + 

-  V  '  ( TUV  ~  Tdv)  ,  (6) 

where  usual  notations  are  adopted.  The  filtered  burn¬ 
ing  rate  Cjv  plus  two  unresolved  flux  terms  appear: 

=  pup-pup,  (7) 

tdv  =  pDpS/ip-pD^fflVip.  (8) 

To  perform  LES,  tu  „  and  tdv  must  be  expressed  from 
quantities  resolved  over  the  coarse  mesh.  Similarly, 
because  of  the  non-linear  character  of  chemistry  4 


wv((p),  and  a  specific  strategy  must  be  developed  to 
express  the  filtered  chemical  source. 

All  these  terms  may  be  computed  from  deconvolu¬ 
tion  operations  providing  approximations  of  unfiltered 
values,  used  to  compute  terms  in  their  primitive  form 
before  they  are  filtered.  In  other  words,  Eq.  (5)  is  ap¬ 
plied  to  get  an  approximation  of  ip(x,  £),  this  approx¬ 
imation  is  used  to  compute  unresolved  terms,  which 
are  subsequently  filtered  using  Eq.  (2). 

For  instance,  the  chemical  source  may  be  estimated 
with  the  mass  fractions  and  temperature  obtained  with 
deconvolution,  to  be  then  explicitly  filtered: 


( x ,  t)  =  U)v  ( C [<£_{x,  t )] ) .  (9) 

In  practice,  <p(x,  t)  is  resolved  by  the  LES  mesh, 
but  the  signal  u>v(C^\y>(x,t)\)  may  not  be  fully  re¬ 
solved  by  this  mesh  and  a  specific  numerical  treatment 
needed. 


Figure  1:  Temperature  profiles.  From  Eq.  (10)  Line:  T{x)\ 
Dash-line:  T(x);  Line  with  circle:  T(x).  From 
solving  Eq.  (15);  Line  with  triangle:  T(x). 
Square:  T(x)  from  deconvolution.  A/5fi  =  8. 


3  Canonical  one-dimensional  fil¬ 
tered  flame  with  Arrhenius  re¬ 
action  rate 

Before  assessing  the  prediction  capabilities  of  approx¬ 
imate  deconvolution  in  a  three-dimensional  turbulent 
flame  with  tabulated  detailed  chemistry,  it  is  verified 
that  simulating  a  planar  and  steady  filtered  flame  pre¬ 
serves  the  flame  major  properties,  as  flame  speed  and 
flame  thickness,  but  also  profiles  filtered  across  the 
flame  for  various  filter  sizes.  A  single-step  Arrhenius 
chemistry  is  considered,  to  test  the  possibility  of  di¬ 
rectly  solving  for  chemistry,  i.e.  avoiding  chemistry 
tabulation  in  future  LES  of  3D  flows. 

The  one -dimensional  canonical  problem  is  cast  in 


terms  of  the  temperature  equation: 

^+V-  (puT)  =  V-(ar(T)VT)+wT,  (10) 
at 

with 

Wr  =  K0(Th  -  T0)pYfYq  exp  ,  dD 

where  K0  =  170  •  1012  s"1,  YF  =  ZSYF>0{  1  -  c), 
Y0  =  (1  -  ZS)Y0,0(  1  -  c)  with  c  =  (T  -  T0)/(Tb  - 
T0),  Zs  =  0.055,  YFo  =  1,  Y0,0  =  0.23  and 
Tac  =  19105  K.  The  Sutherland  law  for  the  diffu¬ 
sion  coefficient  response  versus  temperature  is  intro¬ 
duced:  aT  =  aTo  (1  +  Cs/Ts)/{  1  +  Cs/Tjy/HTjT?) 
with  aFo  =  2.086  •  10-5  m2  •  s-1,  Cs  =  120  and 
Ts  =  291.15  K.  The  fresh  and  burnt  gases  tempera¬ 
tures  are  taken  as:  Ta  =  300  K  and  1),  =  2200  K,  re¬ 
spectively.  The  one-dimensional  unstrained  flame  then 
propagates  at  SF  =  0.47  m  •  s_1;  hence  for  this  pre¬ 
liminary  test,  the  parameters  have  been  chosen  to  pro¬ 
duce  a  laminar  flame  representative  of  simple  hydro¬ 
carbon  burning  in  oxidizer,  under  the  stoichiometric 
condition. 

Approximate  deconvolution/filtering  for  fluxes  and 
sources 

Assuming  constant  pressure  and  molar  mass,  the 
filtered  density  only  depends  on  the  filtered  tempera¬ 
ture  pT  =  JjT  =  p0T0.  The  mass  conservation  reads 
PoSl  =~pu  =  pu,  and  the  velocity  u  =  ( T/T0)SF 
is  also  a  function  of  T.  The  SGS  convective  term  be¬ 
comes  for  this  simple  planar  flame: 

tUt  =  puT-puf ,  (12) 

=  p0SL(T  —  T)  ,  (13) 

=  PoSL  (£^[f]-T)  .  (14) 

Similarly,  the  SGS  diffusive  term  tf>t  is  obtained 
from  deconvolution  and  explicit  filtering,  leading  to 
the  closed  form  of  Eq.  (6)  for  the  temperature: 
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To  check  the  validity  of  approximate  deconvolu¬ 
tion/explicit  filtering  in  this  canonical  case,  the  follow¬ 
ing  procedure  is  applied.  (I):  the  primitive  equation 
(10)  is  solved  with  a  sixth  order  PADE  scheme  (Lele 
1992)  for  spatial  discretization  and  third  order  com¬ 
pact  Runge-Kutta  for  time-stepping.  This  provides  a 


(a)  A/ in  =  3. 


(b)  A/<5fl  =  8. 


Figure  2:  Terms  of  Eq.  (15)  vs  filtered  progress  variable. 

Dash-dot:  Resolved  convection  (i).  Line:  Re¬ 
solved  diffusion  ( ii ).  Circle:  Chemical  source 
(Hi).  Triangle:  SGS  convection  ( iv ).  Dash:  SGS 
molecular  diffusion  (v).  Square:  Total  molecular 
diffusion  (ii)  +  (v).  Dotted:  budget.  Quantities 
are  normalized  by  (SlT0) /Sr. 


reference  laminar  flame  with  a  high-resolution  h  = 
5fi/ 30,  where  5fi  is  the  thermal  flame  thickness.  (II): 
the  modeled  Eq.  (15)  is  solved  with  the  same  numeri¬ 
cal  methods,  but  with  a  mesh  spacing  h  =  A/2,  where 
A  is  the  filter  size.  A  second  order  interpolation  is  then 
added  between  the  coarse  mesh  points  when  comput¬ 
ing  and  filtering  the  chemical  source  from  deconvo- 
luted  temperature.  (Ill):  the  solution  of  (II)  is  com¬ 
pared  against  the  reference  solution  of  (I)  Gaussian  fil¬ 
tered  at  A. 

The  deconvolution  procedure  followed  by  explicit 
filtering  applied  in  Eq.  (15)  exactly  reproduces  the 
Gaussian  filtered  flame,  as  seen  in  Fig.  1  (without  any 
sub-model,  the  flame  cannot  be  converged  over  the 
coarse  mesh).  Tests  were  performed  for  a  range  of  fil¬ 
ters  with  similar  agreement.  In  these  simulations,  the 
flame  speed  computed  from  the  integral  of  the  burning 


Position  (mm) 


Figure  3:  Filtered  temperature  profiles  T(x),  A /Sfi  =  3. 

Bold-line:  Reference  Eq.  (10).  Triangle:  Eq.  (15). 

Square:  Eq.  (17)  with  T(i r)  from  Eq.  (16).  Dash¬ 
line:  Fixed  value  T  =  3. 

rate  is  also  preserved  within  a  few  percent  of  accu¬ 
racy.  Figure  1  also  shows  T(x)  obtained  from  filtering 
of  the  reference  flame  to  illustrate  the  shift  between 
T  and  T,  due  to  density  variation.  In  this  canoni¬ 
cal  one-dimensional  problem,  the  deconvoluted  signal, 
£^1[T](a:)  (square  in  Fig.  1)  perfectly  matches  T(x), 
the  original  one  (line). 

Figure  2  presents  all  the  terms  of  Eq.  (15)  con¬ 
tributing  to  the  temperature  time  derivative.  As  ex¬ 
pected,  the  chemical  source  (circle  in  Fig.  2),  which 
peaks  at  c  =  0.75  in  the  non-filtered  flame,  thus  on 
the  burnt  gas  side,  moves  towards  c  =  0.5  with  mod¬ 
erate  filtering  and  even  on  the  fresh  gas  side  for  larger 
filter  widths  (Fig.  2(b)),  a  shift  that  is  reinforced  here 
by  the  fact  that  the  terms  have  been  divided  by  filtered 
density  to  examine  the  time  derivative  of  temperature 
(Eq.  (15)).  The  SGS  convective  term  (Triangle)  is  pos¬ 
itive  in  the  fresh  gas  and  negative  in  the  burnt  gases, 
leading  to  the  so-called  counter-gradient  contribution 
in  burnt  gases  (Bray  1996;  Veynante,  Trouve,  Bray, 
and  Mantel  1997).  It  is  seen  that  the  SGS  molecu¬ 
lar  diffusion  (Dash-line)  is  of  the  same  order  than  the 
SGS  convective  term  in  this  filtered  laminar  flame  and 
does  not  constitute  a  negligible  contribution  to  the  sum 
of  the  molecular  diffusion  terms  (Square).  The  exact 
contribution  of  this  term  in  a  turbulent  flame  is  still 
to  be  investigated  in  detail.  Increasing  the  filter  size 
decreases  peak  values  of  fluxes  and  sources. 

Approximate  deconvolution/filtering  for  source  and 
corrective  factor  to  resolved  diffusive  budget 

In  a  turbulent  jet  flame,  where  curvature  and  strain 
are  wrinkling  the  flame  front,  the  modeling  specific 
to  one-dimensional  filtered  flames  discussed  above 
may  not  always  be  transposed.  Instead,  a  less  one¬ 
dimensional  oriented  approach  can  be  used  for  unre¬ 
solved  fluxes,  but  still  using  the  deconvolution  for  the 
burning  rate  and  ensuring  that  the  laminar  flame  speed 
is  preserved  whatever  the  mesh  and  the  scalar  signal. 


In  the  thickened  flame  approach  (Colin  et  al. 
2000),  it  was  demonstrated  that  the  correct  propa¬ 
gation  speed  of  a  flame  front  is  secured  applying 
an  amplification  (thickening)  factor  T  to  the  diffu¬ 
sive  budget  expressed  from  mesh  node  values  ( i.e . 
pax^/T  =  FJ>clt(T)X/T)  and  the  same  factor  to  re¬ 
duce  the  source  (i.e.  u>t  =  wt(T) / T),  with  T  >  1. 

As  discussed  in  (Merlin  et  al.  2012),  the  applica¬ 
tion  of  a  corrective  factor  to  the  diffusive  fluxes  may 
be  generalized  to  any  modeled  expression  of  the  burn¬ 
ing  rate  lot,  leading  to  a  dynamically  computed  factor 
T  =  ujt(T)/ujt-  Within  the  present  approximate  de¬ 
convolution/filtering,  the  factor  T  may  be  computed 
according  to  the  explicitly  filtered  source  term,  to  dif¬ 
fer  from  unity  only  in  the  reaction  zone: 


Hx,t)=  IL 

wt(£a  t)\) 


(16) 


The  closed  equation  for  T  then  simply  reads: 


^=-u.VT+iv-  (pTaT{T)\7T'Sj  +  UJT<yC^ 

(17) 

Figure  3  shows  the  temperature  distribution  obtained 
with  this  last  equation.  With  T  computed  dynami¬ 
cally  from  Eq.  (16),  the  response  is  very  close  to  ex¬ 
act  filtering  or  to  the  solution  of  the  full  deconvolution 
(Eq.  (15)),  with  only  a  small  departure  on  the  fresh 
gas  side.  For  comparison,  a  solution  was  obtained  us¬ 
ing  a  fixed  value  of  T  and  lij-  =  /  J7  (TFLES 

model  (Colin  et  al.  2000)).  This  last  option  preserves 
the  flame  speed,  but  cannot  reproduce  carefully  the  fil¬ 
tered  profile  (Fig.  3  dash-line),  which  may  become  a 
limitation  when  applied  to  the  simulation  of  interme¬ 
diate  species  involved  in  more  than  a  single-step  reac¬ 
tion. 

This  preliminary  test  of  approximate  deconvolu¬ 
tion  and  explicit  filtering  suggests  that  the  method  is 
robust  to  simulate  with  precision  the  behavior  of  a  pla¬ 
nar  filtered  flame  front.  However,  this  is  only  a  canon¬ 
ical  problem  and  the  method  is  now  applied  to  the  fil¬ 
tered  chemical  reaction  rate  of  a  fully  turbulent  pre¬ 
mixed  flame,  in  order  to  fully  access  its  prediction  ca¬ 
pabilities. 


4  Turbulent  jet-flame  with  tabu¬ 
lated  detailed  chemistry 

The  premixed  stoichiometric  turbulent  methane  flame 
F3  experimentally  studied  by  Chen  et  al.  (1996)  is 
simulated  (see  Figs.  4  and  5).  The  piloted  Bunsen 
burner  nozzle  diameter  is  D  =  12  mm  and  the  F3 
case  features  a  mean  nozzle  exit  velocity  of  30  m/s, 
with  a  level  of  turbulent  kinetic  energy  of  the  order 
of  3.82  m2/s2.  The  temperature  of  the  stoichiomet¬ 
ric  methane-air  mixture  is  T&  =  2248  K  and  the  pi¬ 
lot  temperature  is  1810  K.  Laser  diagnostics  of  the 
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Figure  4:  Contour  of  iso-Q  criterion  (Q  =  160,  000)  colored 
by  temperature. 


flow  field  are  available,  they  were  obtained  using  two- 
component  and  two-point  laser  Doppler  anenome- 
ter  (LDA).  Scalar  fields  were  also  experimentally  re¬ 
ported  from  2-D  Rayleigh  thermometry  and  line  Ra¬ 
man/Rayleigh  laser-induced  predissociation  fluores¬ 
cence  (LIPF)-OH  techniques. 

The  LES  mesh  is  composed  of  306x194x194 
nodes  over  a  16Dx8Dx8D  computational  domain, 
leading  to  a  resolution  varying  between  170  qm  and 
640  pm.  The  simulations  are  performed  using  the  SiT- 
Com  parallel  flow  solver  already  used  for  testing  SGS 
combustion  closures  (Subramanian  et  al.  2010).  This 
solver  is  based  on  an  explicit  Finite  Volumes  scheme 
for  Cartesian  grids. 

The  Navier-Stokes  equations  are  in  their  fully  com¬ 
pressible  form  together  with  scalars  balance  equa¬ 
tions.  The  convective  terms  are  computed  resort¬ 
ing  to  a  fourth-order  centered  skew-symmetric -like 
scheme  (Ducros  et  al.  2000),  while  the  diffusive  terms 
are  approximated  with  a  fourth-order  centered  scheme. 
Time  is  advanced  explicitly  with  a  third  order  Runge- 
Kutta  method  and  the  3D-NSCBC  boundary  condi¬ 
tions  (Lodato  et  al.  2008)  are  used  at  inlet  and  outlet. 

The  SGS  momentum  fluxes  are  expressed  with  the 
Vreman  model  (Vreman  2004).  To  resolve  eventual 
stiff  gradients  promoted  by  deconvolution,  an  artificial 
dissipation  scheme  was  added  introducing  fourth  order 
dissipative  terms  (Tatsumi  et  al.  1995). 

The  inlet  profiles  are  set  from  the  measurements 
to  which  synthetic  turbulence  is  added  (Klein  et  al. 
2002).  The  total  energy  inlet  condition  accounts  for 
the  non-adiabatic  character  of  the  pilot  flame. 

The  approximate  deconvolution  followed  by  ex¬ 
plicit  filtering  may  be  applied  to  tabulated  chem¬ 
istry.  Considering  a  lookup  table  of  scalars  (chemi- 


cal  sources,  transport  properties,  mass  fractions)  <p  = 
<p(Yc),  with  Yc  the  control  parameters  of  the  tabula¬ 
tion,  the  filtered  scalars  may  be  written: 

W  =  £*$¥>&  [?c]),  (18) 

where  the  relation  (3)  is  used  in  Eq.  (5). 

In  the  present  case,  the  chemical  lookup  table 
is  build  from  laminar  flamelets  computed  with  the 
GRI-mechanism  for  methane-air  combustion  (Smith 
et  al.  1999)  and  the  progress  variable  Yc  is  defined 
as  in  Godel  et  al.  (2008). 

The  equations  solved  for  the  progress  variable  Yc 
and  total  energy  are  similar  to  Eq.  (17),  with  the  chem¬ 
ical  source  obtained  from  the  deconvolution/explicit 
filtering  and  the  factor  T (x,  t)  is  dynamically  deter¬ 
mined  by  the  relation  (16)  applied  to  the  Yc  field. 

In  other  words,  instead  of  explicit  filtering  the  dif¬ 
fusive  term,  the  factor  T (a;,  t )  is  applied  to  the  dif¬ 
fusive  term,  whereas  the  source  term  is  still  obtained 
from  deconvolution  as  in  §  3: 

^  +  V  •  (p5Fc)  =  V  •  (pTDYc  (Yc)  VFC) 

+  (19) 
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Figure  6:  Symbol:  measurements  (Chen  et  al.  1996).  Line: 

LES.  (a)-Left:  Turbulent  kinetic  energy  normal¬ 
ized  by  inlet  value.  (a)-Right:  Time  averaged 
streamwise  velocity  normalized  by  inlet  bulk  ve¬ 
locity. 


Figure  5:  Time  sequence  of  iso-cbc,  filtered  source  of 
progress  variable.  At  =  1.8  ms. 


The  coherent  flow  structure  are  visualized  in  Fig.  4 
with  the  iso-Q  criterion.  The  Kelvin-Helmholtz  large 
scale  vortices  are  rapidly  connected  by  streamwise  in¬ 
stabilities  and  alternate  paring  leading  to  spiraling  vor¬ 
tices,  these  vortices  feature  a  temperature  of  about 
1800  K,  the  coflow  temperature,  from  nozzle  up  to  fur¬ 
ther  downstream.  The  time  sequence  of  iso-burning 
rate  in  Fig.  5  illustrates  the  quite  high  level  of  flame 
wrinkling,  with  even  formation  of  pockets  at  the  flame 
tip.  The  flame  wrinkling  may  be  estimated  in  the  sim¬ 
ulation  as: 


JVc[ 
I  Vc| 


IVgjfl 


(20) 


The  average  value  of  H  computed  over  the  flame  zone 
(0.01  <  c  <  0.99)  at  a  given  instant  in  time  is  (S)  = 
1.1,  with  a  few  local  fluctuations  up  to  3  =  2,  thus 
confirming  the  high-level  of  flame  resolution. 

Statistics  have  been  collected  up  to  convergence 
and  they  are  compared  against  experiments.  The  time 
average  streamwise  velocity  and  kinetic  energy  are 
captured  (Fig.  6),  with  some  deviation  on  the  axis  for 
the  latter.  However,  the  inflow  forcing  has  not  been 
tuned  in  these  simulations,  instead  it  was  kept  fixed 
to  the  profiles  provided  in  the  measurements  at  burner 
exit  plane. 

The  measurements  error  range  reported  by  the 
experimentalists  is  between  8%  and  15%  for  ma¬ 
jor  species  concentration  and  up  to  25%  for  minor 
species  (Chen,  Peters,  Schneemann,  Wruck,  Renz, 


Figure  7:  Mass  fractions  in  %.  Symbol:  measurements. 

Line:  LES.  Left:  Red:  COxlO.  Black:  H2xl00. 
Blue:  OHx75.  Right:  Blue:  O2.  Green:  CH4. 
Red:  H20.  Black:  C02. 

and  Mansour  1996).  Major  species  O2,  CH4,  HoO 
and  C02,  mainly  agree  with  measurements  within  this 
range  (Fig.  7).  The  prediction  of  minor  species  CO, 
H2  and  OH  is  a  sensitive  point  using  flamelet  tabula¬ 
tion  that  does  not  include  strain  rate  as  a  parameter. 
Results  follow  the  experimental  trend  for  these  minor 
species,  with  the  largest  departure  observed  with  H2 
(Fig.  1). 


5  Summary 

A  novel  approach  is  discussed  for  Large  Eddy  Simu¬ 
lation  (LES)  of  premixed  turbulent  combustion,  which 
is  valid  with  high  order  numerics  along  with  grid  res¬ 
olution  providing  quite  resolved  scalar  signals.  It  is 


based  on  an  approximate  and  discrete  deconvolution 
of  density  weighted  scalars  so  that  unclosed  terms  can 
be  estimated  on  the  LES  mesh,  to  be  subsequently  ex¬ 
plicitly  filtered. 

The  method  is  first  tested  for  a  one-dimensional 
laminar  filtered  flame,  using  a  single-step  Arrhenius 
reaction  rate.  Two  options  are  proposed  and  compared 
for  the  fluxes,  a  full  deconvolution,  or,  the  application 
of  a  dynamically  computed  corrective  factor  to  the  dif¬ 
fusive  flux  computed  from  the  resolved  scalar  signal. 

Then,  the  approximate  deconvolution/filtering  of 
the  burning  rate  is  applied  to  three-dimensional  LES 
of  a  turbulent  Bunsen  flame  (Chen  et  al.  1996),  us¬ 
ing  tabulated  detailed  chemistry.  Most  of  the  statisti¬ 
cal  properties  of  both  velocities  and  scalars  are  repro¬ 
duced.  The  modeling  procedure  is  free  from  any  ad¬ 
justable  parameter,  aside  from  numerical  interpolation 
and/or  high-order  damping  of  scalar  gradients,  and  can 
be  readily  applied  to  any  chemical  scheme,  opening 
new  perspective  for  turbulent  combustion  modeling, 
specifically  to  address  intermediate  species  and  pol¬ 
lution  with  LES. 
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